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Abstract 

The acoustical performances of regular arrays of cylindrical ele- 
ments with their axes aligned and parallel to a ground plane have been 
investigated through predictions and laboratory experiments. Semi- 
analytical predictions based on multiple scattering theory and numer- 
ical simulations based on a Boundary Element formulation have been 
made. In an anechoic chamber, arrays of (a) cylindrical acoustically- 
rigid scatterers (PVC pipes) and (b) thin elastic shells have been in- 
stalled with their axis parallel to ground planes consisting either of 
Medium Density Fibreboard (MDF) plate or a sheet of partially reticu- 
lated polyurethane foam. Measurements of Insertion Loss (IL) spectra 
due to the arrays have been made without and with ground planes for 
several receiver heights. The data have been compared with predic- 
tions and numerical simulations. The minima in the excess attenuation 
spectrum due to the ground alone resulting from destructive interfer- 
ence between direct and ground-reflected sound waves, tend to have an 
adverse influence on the band gaps related to a periodic array in the 
free field when these two effects coincide. On the other hand, the pres- 
ence of rigid ground may result in an IL for an array near the ground 
similar to or, in the case of the first Bragg band gap, greater than that 
resulting from a double array, equivalent to the original array plus its 
ground plane mirror image, in the free field. 
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1 Introduction 



Periodic arrangements of acoustic scatterers embedded in a medium with 
different physical properties give rise to band gaps i.e ranges of frequencies 
in which the transmission of acoustic waves is forbidden. If the scatter- 
ers are sohd and the embedding medium is air then these arrays are cahed 
Sonic Crystals (SC). There is interest in the potential use of sonic crystals 
as environmental noise barriers. A semi- analytical approach for predicting 
the transmission properties of sonic crystals has been developed for circular 
scatterer cross-sections and it is based on the superposition of the solution 
for a single scatterer [Tl|2]. However, this scattering approach predicts their 
acoustical performance in the absence of a ground plane. Clearly this will be 
unrealistic if SCs are to be used as noise barriers since a ground will always 
be present. Although the most interesting situation is likely to involve peri- 
odic vertical finite cylinder arrays above a ground plane, this would require 
solution of a 3D problem and hence involve numerical methods and high 
computation resources. Here is considered the more tractable 2D problem 
involving a periodic array of cylinders with their axes parallel to the ground. 

If the ground can be considered to be acoustically-rigid then the mul- 
tiple scattering method can be modified using the method of images to 
construct the reflected acoustic field [3]. For finite impedance ground, it is 
necessary also to satisfy impedance boundary conditions on the ground by, 
for example, using the Weyl-Van der Pol formula [4]. Alternatively, a semi- 
analytical solution has been developed for electromagnetic wave propagation 
that involves an integral representation of the reflected field [5]. Numerical 
approaches can allow for more complex geometries. The Boundary Element 
Method (BEM) based on the integral equation method is the most common 
of these. Specifically it is possible to modify the Green's function [6j so that 
the domain with impedance ground transforms into an unbounded acous- 
tic medium. The result is that the boundary integral equations are only 
considered over the surface of the scatterers. With this approach the com- 
putation time can be relatively low compared to that for the full problem 
with the ground as an additional surface. The method has been widely used 
to predict the performance of the noise barriers in the presence of a finite 
impedance ground [318]. 

In this paper, semi-analytical and numerical methods are used to predict 
the performance of 5 x 3 and 7x3 square lattice arrays consisting of either 
rigid or elastic cylinders with their axes parallel to the ground. The predicted 
performance of these arrays in the presence of rigid or impedance ground is 
compared with their predicted performance in the free field. Insertion Loss 
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Figure 1: Square lattice array above a perfectly reflecting plane, (a) Set of 
vectors used in equation (b) Set of vectors employed in equation ([8]). 

data from experiments carried out in an anechoic chamber are compared 
with predictions. It is shown how the presence of the impedance ground 
affects the IL peaks associated with the so-called band gaps of the sonic 
crystals. 

The analytical and numerical approaches are outlined and some of the 
resulting predictions are discussed in section [2j The experiments are de- 
scribed in section [3l Predictions and data are compared and discussed in 
section |4] before concluding remarks are made in section [5l 

2 Analytical and numerical formulations 

2.1 Multiple scattering 
2.1.1 Rigid scatterers 

Consider a point source and an array of M circular scatterers placed in a 
(positive) half-space characterised by the sound speed in air c = 344 m/s 
and density p = 1.2 kg/m^. Figure [1] illustrates the geometry of the array 
and its image. The position of each scatterer Cm, = 1..M, is given by 
the vector Rm- The position of the image of scatterer Cm is defined by the 
radius vector i?^. The scatterers are considered to be arranged in a square 
lattice which is defined by the lattice constant L. However the methods 
described subsequently can be applied to any other lattice configuration. 
The solution of the appropriate scattering problem satisfies the Helmholtz 
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equation in the half-space that is written in polar coordinates (r, 9) as 

Ap{r) + k^p{r) = 0, (1) 
I d f d\ I 

where A = -— r— H — jttt;^, r = r(cos 6, sin 9) is the radius vector, p is 
r or \ or J afc*^ 

acoustic displacement potential, A; = w/c, w is angular frequency. Equation 

([T]) is solved in conjunction with radiation conditions 

— ikp = o (j"'^^"^^ , as r — )• oo, (2) 

and with the Neumann condition imposed on the boundary of acoustic half- 
space (i.e. rigid ground) and on the surface of the scatterers (this condition 
has to be replaced by continuity conditions if scatterer is an elastic shell 
[ini) that is 

Using the multiple scattering technique [TJ [S] and the method of images [5] 
the general solution of the formulated problem can be written as [11] 

p{r)=Po{r)+Ps{r), (4) 

whereby contributions from the point source and its image are collected in 
Po i-e. 

Po(.r) = po,d{r) + po,r{r), (5a) 
PoA^) = Hl^'\kro), (5b) 
PoAr) = Hi,'\kr'o), (5c) 

whereas scattered direct and reflected acoustic fields are described by 

Ps{r) =PsA'^) +Ps,r{r), (6a) 

M +00 

PsAr) = Yl E ^™^r^P(fcr„)e-^-, (6b) 

m=l n=— cxD 

M +0O 

PsAr)=Yl E A^Z:^Hi'\kr'Je-^-'-. (6c) 

m=l n=—oo 

The vector vq = ro(cos sin^o) connects the point source and the receiver 
point (i.e. point P in Figure [T]J a)). The vector = rm{cos9m,sm9m) 



4 



connects the centre of scatterer Cm and the receiver. The prime (') is used 
to indicate the geometrical parameters for the image source and the image 
scatterers placed in the negative half-space. A!^, n £ 7j, m = 1..M are 
unknown coefficients. 

The factors describe the type of conditions imposed on the surface of the 
scatterers and in case of rigid cylinders they can be expressed as [1] 

rym _ 9rJn{ka.m) 

~7 — (Ty: — r- ^'^ 

where am is the radius of scatterer Cm and dr is the derivative with respect 
to polar coordinate r. 

The solution for the unbounded acoustic space can be retrieved from 
equation (j3|) by putting to zero in equations (f5a|) and ([6a|) all terms related 
to the constructed images that are po,r and ps^r respectively. One can also 
deduce from equation dH) and vector definitions in Figure [l]||a) that for the 
source and receiver both on the ground the acoustic pressure in a half-space 
is double the pressure in the unbounded acoustic space. 

Applying the addition theorem [131ll4j. described in Appendix A, to the 
solution dH) and substituting it to the boundary condition ([3]) , the algebraic 
system of equations can be derived to find the unknown coefficients A^. 
This system is given by 

a::+ (8) 

oo ( M M 
g=— oo I p=l,p^m p=l 

where vector Rom = ^Om(cos aom, sinaom) defines the position of scatterer 
Cm with respect to point source and vector Rmp = Rmp{cos amp, sin amp) 
defines the position of scatterer Cp with respect to scatterer Cm- Again the 
system of equations ([8]) can be transformed to that for the case of unbounded 
acoustic space by eliminating all terms dependent on the geometrical param- 
eters of image source and scatterers. To be solved numerically the infinite 
system of equations (l8|) is truncated to the finite number of M{2N + 1) 
equations. If 5 < < 7, the numerical solution is accurate up to four 
significant figures [1]. It is also noted that for the considered configurations 
and frequency range the computation time required to solve system ([8]) on 
an Intel Core 2 Duo processor based PC is between 60 and 180 s. This is 
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true for the codes executed in the GNU Compiler Collection (GCC) as well 
as for the scripts using a commercial software MATLAB. 




Figure 2: (a) Predicted insertion loss spectra with source and receiver coor- 
dinates of (0,0) and (10,0) respectively (i) for a 5 x 3 square array of rigid 
cylinders with Um = 0.1 m above acoustically-rigid plane at y = 0, (ii) for 
the same array in the free field and (iii) for the array plus its mirror image 
(a 10 X 3 array) in the free field, (b) Diagram showing that the 10 x 3 array 
in the free field consists of the original 5x3 array plus its mirror image in 
the ground plane, (c) Diagram of the 5x3 array over rigid ground. 

In Figure O the predicted insertion loss spectrum due to a 5 x 3 array 
of rigid scatterers over an acoustically-rigid ground is compared with those 
obtained (a) for the same array in free field conditions and (b) for the original 
array plus its mirror image array i.e. a 10 x 3 array in free field conditions. In 
all configurations the nearest part of an array from the source is at = 1.5 
m. Also note that throughout this paper the insertion loss is calculated as 



IL = 20 log 
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Po 

p 



(9) 



In the free field the cylinder locations in the lower half of the 10 x 3 array are 
defined by the coordinates of the image cylinders in the half-space problem. 
The distance to the ground, Hy = 0.15 m, of the centers of the lowest 
cylinders in the array is half of the lattice constant L = 0.3 m so that they 
are separated from the centers of the cylinders of the image array nearest 
the ground plane by the lattice constant. This means that for the geometry 
considered the array and its mirror image in the ground plane effectively 
form a complete regular array of twice the size. With the source and receiver 
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Figure 3: Predicted insertion loss spectra with source coordinates (0,0.235) 
m for 7x3 array of rigid cylinders with am = 0.0275 (i) over acoustically- 
rigid ground with receiver coordinates (1.203,0) m (solid black line) (ii) 
over acoustically-rigid ground with receiver coordinates (1.203,0.235) (bro- 
ken line) and (iii) in free field (dash-dot line) with receiver coordinates 
(1.203,0). 



on the ground, the predicted insertion loss spectrum of the 5x3 array in 
the presence of the rigid ground is the same as that predicted for an array of 
double the size (10 x 3) in the free field. It is also observed that the insertion 
loss of a 5 X 3 array in a half-space is predicted to be higher near 573 Hz 
(the first Bragg band gap) than that for the same size of the array in the 
unbounded acoustic space. 

Figure [3] compares predicted insertion loss spectra for 7x3 array of 
rigid cylinders with its counterpart in the free acoustic field. As before the 
cylinders are arranged in a square lattice with L = 0.069 m. The nearest 
part of the array from the source is at = 0.755 m and the distance of the 
array to the ground is Hy = 0.0345 m. The predicted effect of raising the 
receiver is clearly detrimental to insertion loss at frequencies corresponding 
to the (rigid) ground effect dip. It is also observed that performance of the 
array over the rigid ground with receiver on the ground is predicted to be 
improved between 2000 Hz and 3000 Hz compared to that in the unbounded 
acoustic space. 



2.1.2 Elastic shell scatterers 

A multiple scattering analysis can be carried out to predict the insertion 
loss spectrum due to an array of elastic shells with their axes parallel to a 
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Figure 4: Predicted insertion loss spectra for the source-receiver-array 
(5 X 3) geometry specified for Figure (H^c) in the presence of acoustically- 
rigid ground at y = 0, with rigid cylinders (solid line), elastic shells 
{ps = 1650kg/m3, E = 1.75 MPa, u = 0.4998, C2 = 23 m/s, a„ = 0.1) 
with wall thickness 2h = 0.001 m (broken line) and elastic shells with wall 
thickness 2h = 0.002 m (dot-dash line). 



rigid ground. The identical elastic shells are characterised by their density 
Ps, Young's modulus E, Poisson's ratio i^, shear velocity C2, half-thickness 
h and the mid-surface radius S = — h. For certain ranges of values 
of these parameters, the first elastic shell resonance (i.e. the axisymmetric 
resonance) can be observed below the first Bragg band gap associated with 
the lattice constant of the array in the unbounded acoustic space. This 
results in additional positive insertion loss peaks [10]. 

The asymptotic theory of thin elastic shells [15] has been used \10\ [TT] 
to derive the factors 

^ drJn{kS) 



drHi^\kS) + iUn 



(10) 



where 



Un = - . " (11) 



K TiSh (l + n2 - /t|52) drJn{kS) ■ 



e = pc/ {psC2) is the relative impedance, k = c/c2 and = uj ^JpJ^l^^^J^Yi^ ■ 
If the relative impedance (e) tends to zero, then C/„ becomes negligible and 
the form of in (fTO|) reduces to that in d?]). 

Figure [4] compares the predicted insertion loss spectra of the array of 
elastic shells in the acoustic half-space with that of the array of rigid shells. 
These results are similar to those in the unbounded acoustic space with array 
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plus its mirror image array i.e. doubled in size. Additional insertion loss 
peaks due to axisymmetric resonances of the elastic shells are observed below 
the peak related to the Bragg band gap. The frequency of the axisymmetric 
resonance reduces with the increased shell thickness. 

FigureEJcompares insertion loss spectra for 7x3 rigid cylinder and elastic 
shell arrays above acoustically-rigid ground with different receiver heights. 
Scatterers are arranged in square lattice with L = 0.069 m. The position 
of the arrays with respect to the source and ground plane is identical to 
that described in Figure El It is shown that the existence of the predicted 
effect due to axisymmetric resonances of the elastic shells (see Figure H|) 
is dependent on the geometrical parameters of the problem such as receiver 
coordinates. For the receiver heights 0.117 m and 0.235 m considered in Fig- 
ures ISJa) and (b) respectively, the additional peak due to the axisymmetric 
resonance of the shell appears around 1000 Hz. On the other hand this peak 
does not exist when receiver is raised to the 0.352 m (see Figure [5]^c)) since 
there is a destructive interference in the ground effect at this frequency. 

2.2 Calculations based on the boundary integral equation 

To investigate the influence of finite impedance of a ground plane on the 
insertion loss due to an array of regularly spaced cylinders parallel to the 
impedance surface the solution has been sought to an appropriate boundary 
integral equation. 

The Laplacian in equation ([T|) is rewritten in terms of (x, y) coordinates 



The boundary condition imposed on the ground surface is written as 



where j3 is admittance of the homogeneous impedance plane [6, eq. (1.2.11)]. 

Then, applying relations ([I]),©, ()12p and condition of rigid scatterer 
surface dp/dr = to the Green's theorem [16] the integral equation for p{r) 
can be derived in the following form |1H [12] 



using A = d"^ /dx^ + d"^ /dy"^. 



— ikPp = 0, 



(12) 




(13) 



where 




m 



(14) 
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Figure 5: Predicted insertion loss spectra with source coordinates (0,0.235) 
m for a 7 X 3 array of (i) rigid cylinders of radius = 0.0275 m (solid 
black line) and (ii) elastic shells of radius = 0.0275 m and thickness 
2h = 0.00025 m (broken line) over acoustically-rigid ground with receiver 
coordinates (a) (1.203,0.117) m, (b) (1.203,0.235) m and (c) (1.203,0.352) 
m. The elastic shell is made of latex with material parameters specified in 
Figure HI 
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Figure 6: Predicted insertion loss spectra with source and receiver coordi- 
nates (0,0) m and (10,0.45) m respectively for a square 5x3 array of rigid 
scatterers of radius 0.1 m with lattice constant L = 0.3 m, Hx = 1.5 m and 
Hy = 0.15 m. The array of scatterers is placed above (a) an acoustically- 
rigid ground (scattering theory - solid line; boundary element calculation - 
broken line) and (b) an impedance ground (the solid, broken and dash-dot 
lines correspond to the different values of effective flow resistivity in the 
key). 

with r = {x,y), n(rs) is the unit vector normal to the scatterer surface 
and directed outward, and dCm is the surface of scatterer Cm- G^(ro,r'), 
which is the solution for a half space above an impedance plane, is given by 
equations (2.1.2), (2.1.20), (2.1.21), (2.1.44) and (3.6.21) in [6] and is not 
repeated here. Note, that in relation (jl4|) the corner points of an obstacle 
are not defined due to the circular shape of the scatterers. 

Figure [61(a) demonstrates that the boundary integral formulation yields 
results close to those obtained using multiple scattering theory for an array 
of horizontal cylinders above acoustically-rigid ground. In case of the rigid 
ground the computation time is comparable with that of the semi-analytical 
method ([8]). The difference in the predictions at higher frequencies can be 
reduced by finer discretization of the surface of the scatterers. This however 
increases the computation time. 

Figure EJb) shows predictions obtained using the boundary integral for- 
mulation for three values of ground impedance based on a one parameter 
(effective flow resistivity) impedance model |17l HI]. It is noted that the 
computation time required to perform the numerical calculations for the 
impedance ground although less than an hour is substantially bigger than 
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that for the rigid ground. The predicted insertion loss spectrum for the 
lowest value of effective flow resistivity (20 kPas/m-^ corresponding to a 
mineral wool) shows more or less complete elimination the band-gap effect 
whereas the predicted insertion loss spectra for the higher flow resistivities 
(168 kPas/m^ and 250 kPas/m^ corresponding to hay and grassland respec- 
tively) indicate that in the presence of a relatively acoustically-rigid surface 
the XL spectrum due to the 5x3 array is predicted to include maxima in 
the frequency intervals corresponding to the array band-gaps. 

An alternative approach to BEM is that based on the Weyl-Van der 
Pol formula modeling locally reacting ground [19J. This approach has been 
employed for a single scatterer above an impedance plane [20^ I21j . Com- 
pared to BEM the use of the Weyl-Van der Pol formula has the advantage 
of reduced computation time comparable with that of the semi-analytical 
approach for the scatterers over the rigid ground. However its application 
to an array of scatterers is heuristic and can only be used within a limited 
range of source-array and array-receiver distances. An example in section H] 
shows that the results deteriorate with increasing receiver height. 

3 Laboratory experiments 

Measurements of the insertion loss spectra due to arrays of regularly spaced 
parallel rigid cylinders and elastic shells without and with ground planes 
have been carried out in an anechoic chamber. Rigid cylinders consisted of 
2 m long PVC pipes with outer diameter 0.055 m. 2 m long elastic shells 
were made from 0.25 mm thick sheets of Latex by overlapping the edges 
and gluing them together. The sound source was a Bruel & Kjaer point 
source loudspeaker controlled by a Maximum-Length Sequence System An- 
alyzer (MLSSA) system enabling determination of impulse responses in the 
presence of noise. A Bruel & Kjaer 1/2 inch microphone was used as the re- 
ceiver. Figures [7] (a) , (b) and (c) show example measurement arrangements. 
Supports for the 2 m long cylinders were provided by holed MDF boards at 
the top and base of each array. To maintain their shape and vertical orien- 
tation, the latex cylinders were slightly inflated above atmospheric pressure 
through a common pipe connecting to a small pump. 

An MDF board large enough to avoid the diffraction at the edges was 
used as a rigid surface. The rigid cylinders could be arranged horizontally 
above a horizontal MDF board (Figure [7|^b)). However since the latex cylin- 
ders had to be arranged vertically to preserve their shape, the MDF board 
was also supported vertically (Figure [T^c)). In both cases the cylinder axes 
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(a) (b) (c) 

Figure 7: (a) Schematic of the experimental configuration showing the source 
location, the sonic crystal array and three receiver locations, (b) Photograph 
of experimental arrangement with rigid cylinder array above an impedance 
(polyurethane foam layer) ground, (c) Photograph of experimental arrange- 
ment with Latex shell array near to rigid (MDF board) ground 



were parallel to the board. For simplicity all distances to the MDF board 
in both setups are referred to as heights. As shown in Figure [Tl^a), the loud- 
speaker point source was positioned di = 0.755 m from the array of rigid 
cylinders and di = 0.35 m from the array of Latex shells at the height of 
the horizontal mid-plane of the array (hs = 0.23 m above the ground). The 
height of the receiver microphone was hri = 0.117 m, hr2 = 0.235 m or 
hr3 = 0.352 m and it was placed in a vertical plane ds = 0.257 m from the 
back of the array. The receiver heights were chosen to be below, at, and 
above, the horizontal mid-plane of the array. In all cases, the distance be- 
tween the microphone and the cylinder array has been considered the same. 
The difference between the sound levels recorded in the X direction (0°) at 
the same point behind the array with and without the ground was measured 



4 Comparisons between data and predictions 

Figure [8] compares measured and predicted insertion loss spectra for a 7 x 3 
rigid cylinder array over rigid ground for three receiver heights. The predic- 
tions assume the source-array-receiver geometries used in the experiments 
described in section [3l Up to 1500 Hz the predictions and data are in close 
agreement. Above 1500 Hz there are some discrepancies which may be due 
to unwanted reflections and departures from the assumed ideal geometry. 
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Figure 8: Measured (solid line) and predicted (broken line) insertion loss 
spectra for a square 7x3 array of rigid cylinders of diameter 0.055 m over 
acoustically-rigid ground with source coordinates (0,0.235) m and receiver 
coordinates (a) (1.203,0.117) m, (b) (1.203,0.235) m and (c) (1.203,0.352) 
m. 

Both data and predictions in Figures [8{a) and [8]^c) show IL maxima near 
2500 Hz which are associated with the Bragg band gaps expected in the 
unbounded domain. Both data and predictions for the elevated receiver 
heights (0.235 m and 0.352 m) show the adverse influences of destructive 
interference associated with the (rigid) ground effect on the IL spectra near 
2000 Hz and 1250 Hz. Indeed in Figure [8|^b), for the receiver at 0.235 m 
height, it is clear that the destructive interference in the ground effect near 
2500 Hz is rather dominant. 

The measured and predicted performances of 7 x 3 array of Latex shells 
over the rigid ground in Figure [9] are similar to that described in section 
12.1.21 In particular. Figure [9]|^c) shows that when the first ground effect dip 
is in the vicinity of the axisymmetric resonance of the shell (900 Hz) the 
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Figure 9: Measured (solid line) and predicted (broken line) insertion loss 
spectra due to a square 7x3 array of Latex shells of diameter 0.055 m, 
thickness 0.00025 m and material parameters specified for Figure U] over 
acoustically-rigid ground. The source is at coordinates (0,0.235) m and the 
receiver coordinates are (a) (1.203,0.117) m, (b) (1.203,0.235) m and (c) 
(1.203,0.352) m. 

corresponding positive IL peak is no longer present. 

Figure [10] compares the measured and predicted insertion loss spectra 
for 7x3 rigid cylinder arrays over finite impedance ground using the source- 
array-receiver geometries described in section [3l To obtain the predictions 
in Figure [10] the properties of impedance of the hard-backed foam layer 
have been deduced from best fits short range measurements of complex 
excess attenuation [23j. As a result the finite impedance (open cell foam 
layer) surface is represented by a two parameter impedance model with cje = 
4kPas/m^, = 105 m~^. There are discrepancies between predictions and 
data over the whole frequency range but the predictions follow the general 
trends in the data. Compared to the results for the rigid ground plotted in 
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Figure [8] the IL minima associated with the ground effect are shifted towards 
lower frequencies. Both measurements and predictions show that as a result 
of this shift in the ground effect the IL maxima related to the Bragg band 
gap can be observed for all three positions of the receiver. 




Frequency — Hz 
(c) 27r 



Figure 10: Measured (solid lines) and predicted (broken lines) insertion loss 
spectra for a square 7x3 array of rigid cylinders of diameter 0.055 m over 
finite impedance ground with source coordinates (0,0.235) m and receiver 
coordinates (a) (1.203,0.117) m, (b) (1.203,0.235) m and (c) (1.203,0.352) 
m. 

The Weyl-Van der Pol formula [19] can also be used to predict the field 
due to a point source above an impedance plane by using the multiple scat- 
tering technique and method of images. For a line source over an impedance 
ground the acoustic wave field is approximated by [21J 

Po{r) = H^^\kr^) + QQH^^\kr'^), (15) 

where Qq is the spherical wave reflection coefficient described below by equa- 
tion ([IT]). 
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Figure 11: Measured (solid line), predicted with BEM (broken line) and pre- 
dicted with the Weyl-Van der Pol forniTila (dash-dot line with open circles) 
insertion loss spectra for a square 7x3 array of rigid cylinders of diame- 
ter 0.055 m over finite impedance ground with source coordinates (0,0.235) 
and receiver coordinates (a) (1.203,0.117) m, (b) (1.203,0.235) m and (c) 
(1.203,0.352) m. 



By the analogy with the source over the ground the scattered wave field 
for the array of circular scatterers can be written as 

M +00 
m=l n=— oo 

The spherical wave reflection coefficient Qm, rn = 0..M, is given by 

Qm = Vm + il- Vm)F{Wm) (17) 
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where 



Vm = 



COS am — P 



(18a) 



cos + /? 




(18b) 



m 



) = 1 + i^/nwm exp (--w^) erfc(-mm) 



(18c) 



within which Om is the angle of incidence defined by either position of the 
source or centre of the scatterer [21] . 

In Figure [TT] predictions based on (a) BEM and (b) the semi-analytical 
approach described by equations ([H) and p^ - (fT8]) are compared with the 
measured insertion loss for 7x3 array of rigid scatterers over the impedance 
ground with the parameters identical to those used for Figure [TOl Figure 
\TT\ a) shows that semi-analytical results using p^ - (fT8]) are in good agree- 
ment with the data if the receiver is close to the ground. However, when the 
receiver is at heights of 0.235 m and 0.352 m, the semi-analytical approach 
predicts a peak and dip in the IL above 4000 Hz and a dip near 2500 Hz 
respectively that are not observed in the data or in the BEM predictions 
(see Figures flTT b) and [TlTc)). It may be concluded that the accuracy of 
the heuristic semi- analytical solution (equations (fT5]l - (fT8]) ) decreases as the 
receiver height increases. 

5 Concluding remarks 

Semi-analytical and numerical models have been derived for predicting mul- 
tiple scattering effects of a finite arrays of cylinders parallel to rigid and 
impedance ground respectively. The numerical technique (BEM) has been 
validated against the semi-analytical multiple scattering approach for rigid 
cylinders above rigid ground. Results of both methods have been compared 
with data. It has been shown that performance of an array in a half-space is 
similar to that of the doubled array (i.e. an array composed of the original 
array plus an array corresponding to its mirror image in a rigid plane) in the 
unbounded acoustical space subject to conditions reported in the discussion 
of Figure [21 Depending on the source- array-receiver geometry the presence 
of a rigid ground can result in destruction of the positive IL peak associated 
with the first Bragg band gap by the first destructive interference minimum 
in the ground effect. However introduction of the impedance ground results 
in the shift of ground effect minima to lower frequencies so that the Bragg 
band gap is maintained. The numerical BEM technique for predicting the 
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IL spectra due to finite cylinder arrays over impedance ground has been 
compared with an alternative semi-analytical approach based on the Weyl- 
Van der Pol formula. The results show that heuristic approximation of the 
influence of the impedance ground in the semi-analytical approach becomes 
worse as receiver height is increased. 
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Appendix A: Graf's addition theorem 




Figure 12: Geometry for Graf's addition theorem. 

In this section Graf's addition theorem is modified so that it can be 
applied to the solution of the reflected scattered field ps,r in equation ([4]). 
First the addition theorem is stated for the solution of the direct scattered 
field ps,d, yielding 

oo 
q= — oo 

for rm < Rmp, rn £ Z. The outlined form of the additional theorem is based 
on the configuration shown in Figure [T2l To adapt theorem (jl9p to solution 
Ps^r the index n has to be replaced by its negative counterpart n = —n. 
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Using the relation H_l{z) = e^^'^Hn {z) the addition theorem is written as 

oo 

g=— oo 

The latter can be used in equation @ to transform image solution to that 
defined by the variables of the real scatterer. 
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